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Abstract 

Cosmological consequences of a string-motivated dark energy scenario featuring a scalar field coupled to the Gauss- 
Bonnet invariant are investigated. We study the evolution of the universe in such a model, identifying its key properties. 
The evolution of the homogeneous background and cosmological perturbations, both at large and small scales, are 
calculated. The impact of the coupling on galaxy distributions and the cosmic microwave background is examined. 
We find the coupling provides a mechanism to viably onset the late acceleration, to alleviate the coincidence problem, 
and furthermore to effectively cross the phantom divide at the present while avoiding a Big Rip in the future. We 
show the model could explain the present cosmological observations, and discuss how various astrophysical and 
cosmological data, from the Solar system, supernovae la, cosmic microwave background radiation and large scale 
structure constrain it. 

Key words: Astrophysics and Cosmology 
PACS: 98.80.-k, 98.80.Jk 



Einstein General Relativity together with ordi- 
nary matter, described by the standard model of 
particle physics, cannot fully explain the observa- 
tional data from Supernovae type la (SNela) [1], 
the matter power spectrum of large scale struc- 
ture(LSS) [2] and the anisotropy spectrum of the 
Cosmic Microwave Background Radiation(CMBR) 
[3]. One needs to introduce two exotic components 
into the matter-energy budget of the Universe. Dark 
matter, a fluid with zero or very small pressure, 
corresponds to about 25% of the universe energy 
budget. Dark energy, with its negative pressure, 
dominates the Universe density and is responsible 
for its present acceleration [4] . 
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The observed acceleration and the undisclosed na- 
ture of these two exotic components strongly mo- 
tivates the extension of both the standard model 
of particle and gravitational physics sectors [5,6,7]. 
This is also what quantum gravity seems to require. 
Typically the low-energy limit of string theory fea- 
tures scalar fields and their couplings to various cur- 
vature terms. Interestingly, there is a unique com- 
bination of the curvature squared terms, the Gauss- 
Bonnet (GB) invariant 

^GB = R^^^'^ Rfiupa — ^R'^'^Rfiu + R^, (1) 

that is both ghost-free in Minkowski backgrounds 
and leads to second order order field equations. All 
versions of string theory in 10 dimensions (except 
Type II) include this term as the leading order a' 
correction [8,9]. Specifically, these couplings can be 
shown to appear at the tree level or the one-loop 
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level (depending on whether one considers a dilaton 
or an average volume modulus) of string effective 
action when going from the string frame to the Ein- 
stein frame [10,11]. The effective action could then 
be written as 



k"^^, II = 8K'^(l)Hf'{(f)), and e=^.We consider 
a canonical scalar field with 7=1 and adopt an ex- 
ponential form for the potential V{(f>) = Vqc"'^'^'^/^ 
and for the coupling / = /ge"'*'^/^. The nonpertur- 
bative effects from gaugino condensation or instan- 
tons can result in an exponential potential [19]. An 
exponential field-dependence can approximate the 
coupling ensuing from, for instance heterotic com- 
pactification [10,11]. 

The background has altogether six fixed points 
when considered dynamical system. However, 
only two of them are now relevant (see [20] for de- 
tails). The first is the standard exponential tracking 
solution, where the scalar field mimics exactly the 
background equation of state Wm- This fixed point 
is a stable spiral when A > a, \/6(l -|- Wm)- It is, 
however, a saddle point when A, a < \/6(l + w„i) 
or a > A, •\/6(l + Wm)- The later happens when the 
coupling becomes significant at late times. Then the 
field will be passed from the scaling solution to a 
potential-dominated solution (see [20]). 

Since for the scaling solution we have H"^ r-^ p ^ 
^-3(i+tum)^ the last term in the Fricdmann equation 
scales like pf = 8H^f'{(p)<p ~ a-3(i-H«.„)(2-a/A)_ 
This follows from the tracking behaviour of the 
scalar field; since (f>' = 3(1 + w„i)/\, wc have that 
4> = (po + 3(1 -|- Wm) log(a)/A, and hence f'{(j)) ~ 
^3a;(i-i-w,n)/A Tii^g fln^ tj^^t the effective energy 
density due to the presence of the GB term, pf, di- 
lutes slower than the energy density due to matter, 
Pm, if and only if a > A. 

A typical evolution for the model is shown in FIG. 
1. Note that > 1 may occur, since the (effective) 
GB energy density 0/ can be negative , and that 
Weff < — 1 is possible here as well. 

To investigate the cosmological effects of the GB 
coupling in more detail, we also consider the linear 
perturbations. In the synchronous gauge [21], one 
can parameterize the metric perturbations as 



-I 



R 
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where K = (87rG)-i/2. The scalar field La grangian 
is = —^(V(/))^ — y ((/)), where 7 is a constant. The 
function f{(j)) = a — 5£,{(f)): the coupling a may be 
related to string coupling via a ^ l/ff^. The nu- 
merical coefficient 5 typically depends on the mass- 
less spectrum of every particular model [11]. 

In this work we explore the cosmological conse- 
quences and viability of the action (2), with spe- 
cific choices for the potential and the coupling, both 
taken to be single exponential terms. Such mod- 
els have been recently studied in [12] within the 
context of the dark energy problem, and their pos- 
sible background evolution has been investigated 
[13,14,15,16,17,18]. We wih show that even the sim- 
ple and well-motivated exponential parameteriza- 
tion can naturally exhibit a viable transition to ac- 
celeration and a transient phantom expansion. Be- 
sides presenting such possibility in the model and 
quantitatively testing it with various observational 
data, we also consider the phenomenology of GB 
dark energy models in general at cosmological, as- 
trophysical and Solar system scales during the whole 
cosmological evolution. In particular, we investigate 
how the coupling affects the CMBR anisotropies and 
LSS. 

We consider a flat, homogeneous and isotropic 
background universe with scale factor a{t). Deriva- 
tives with respect to the cosmic time t arc denoted 
by an overdot, and a prime means derivative with 
respect to the e-folding time log (a) unless other 
variable is explicitly specified. Action (2) yields the 
Fricdmann equation 

\H^ = U-' + V{<t>)+Pm + 2AH^rm, (3) 



and the Klein- Gordon equation 

7(<?i + 2,H^) + y (0) + S\cl>)Rls = 0, 



(4) 



where H = d/ais the Hubble rate, p„i represents the 
matter component, and the GB invariant is Rqq = 
2AH'^{H + H'^). It is convenient to define the di- 



-df + a'^{t){dij + hij)dx'dx^. 



Conventionally the scalar modes are then defined 
in the Fourier space by the decomposition into the 
trace (/i) and the traceless part {rf) of hij , 



/iy(x,T) = / d^kf^^^'lkkK^^T) 



mensionless variables Q.n 



6{kikj - ^Sij)r]{k, r) 



where k = /c/c. One can then write the energy con- 
straint (perturbed version of the Fricdmann equa- 
tion) as an evolution equation for the metric poten- 
tial h, 
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Fig. 1. Background evolution when A = 4 and a = 20. In 
general, the evolution in the model consists of 1) the scaling 
attractor and 2) the potential-dominated (de Sitter) solution. 
Existence of 1) requires A > \/6, and the transition to 2) then 
occurs if o > A. Here we plot fractional energy densities for 
matter, (dash-dotted line), scalar field = [x^ +y)/Z 
(dashed line) and the GB term, Qf = fj. (dotted line). Solid 
line is the total equation of state ui^^j = — 2e/3 — 1. 
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where the density fluctuation is as usual Sp = 
PmSm + 4>^4' + V{<j))S(f>. Note there appears both 
new source terms due to the fluctuations in the 
coupling / as well as a modulating prefactor chang- 
ing the response to the standard source termfjJJ. 
Similar modifications are present in the momentum 
constraint equation governing the evolution of the 
other metric potential 77, 



V = 



AttG 



a{p + p)e - 8H{H6'f - H^Sf) 



(6) 



where the velocity perturbation is as usual 
a{p + p)9 — a{pm + Pm)Om + k(j>S(f>. The Klein- 
Gordon equation for the scalar field fluctuation, 
being the perturbed version of Eq.(4), is also non- 
trivially modified [22] . Since still minimally coupled 
to gravity, all other matter obey the usual continu- 
ity equations [20,23]. 

We have numerically integrated the fully per- 
turbed equations and computed the full matter 
power and CMBR spectra for the example model 
presented here. In FIG. 2 we show how the cosmo- 






^ The expression implies divergence when ^ = 2/3 (unless 
it would happen that the square bracket term vanishes for 
all fc-modes at that point). However, in the models consid- 
ered in this letter we always have the relative Gauss-Bonnet 
contribution /i < 2/3. Then also Eq.(6) is well-behaved. 
^ We used a modified version of the CAMB code [31]. 



Fig. 2. Top two figures: The effect of the potential slope on 
the CMB and matter power spectra. Here a = 20. Dotted 
lines are for A = 4.5, dashed line for A = 6.0, and dash- 
dotted for A = 8.0. Bottom two figures: The effect of the 
coupling slope on the CMB and matter power spectra. Here 
A = 6.0. Dotted lines are for a = 10, dashed line for a = 20, 
and dash-dotted for a = 30. The solid line is ACDM model. 

= 0.4 for all fi gures. The CMB and matter power spec- 
tra error bars are from the WMAP data [3] and SDSS[2] 
respectively. 



logical predictions are changed when the slope of 
the potential or of the coupling are varied. The main 
imprint from different potential slopes A seems to 
be in the normalization. For low values of A, there 
is significant contribution of the scalar field during 
the matter dominated era. This slows down the 
rate of growth of matter inhomogeneities. Hence 
the fact that there is less structure nowadays than 
for larger A is not a consequence of the GB modi- 
fication, but rather an effect of the presence of the 
scalar field in the earlier scaling era. Finally, in the 
bottom of FIG. 2 we see that the strength of the 
coupling a might be more difficult to deduce from 
these data. With steep coupling slopes, the scalar 
field domination takes place more rapidly and with 
more negative w^ff, which can somewhat amplify 
the ISW effect. The contrary happens for smaller a. 

We now consider the constraints arising from as- 
trophysical and cosmological observations. We cal- 
culate the SNela luminosity-distance relation. To 
compare with data, we use the " Gold" sample of 157 
SNela from Rcf. [1] and marginalize over the Hub- 
ble constant Hq. We also compute the CMBR shift 
parameter [24] TZ and apply the latest constraints 
[25]. The combined constraints arising from all these 
data is shown in FIG. 3. The SNela data alone is con- 
sistent with a wide range of matter densities flm, 
but when combined with the CMBR parameter TZ 
one is restricted to rather high H.^ ^ 0.4. Note that 
we restrict here only to cosmologies where a scal- 
ing matter era is followed by the acceleration era. 
The existence of the scaling attractor requires A > 
•\/6, otherwise the evolution is different and depends 
on the initial conditions for the field. Therefore we 
leave models with lower potential slopes, a < -s/G, 
out from this study, but their phenomenology could 
be interesting to investigate elsewherJ^. 

Time variation of the effective gravitational con- 
stant is tightly constrained by observations 
within the Solar system and laboratories, and indi- 
cate that |G',/G,| < 0.01 [26]. To derive the varia- 
tion of this constant for the coupled GB gravity, we 
follow the approach of Ref. [27] where cosmologi- 
cal perturbation equations were considered at their 
Newtonian limit. The Poisson equation was derived 
and the effective strength of gravitational coupling 
was read from the resulting expression, which re- 



^ Such potentials are in accordance with various compacti- 
fication schemes and particle physics models. One also notes 
that then o > A is not required for the transition to accel- 
eration. 
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Fig. 3. The 68, 90 and 99 percent confidence limits for the 
model in the Qrn - A plane when a is marginalized over in 
the range 1.5A < a < lOA. Dotted lines are constraints from 
the SNela data and the solid lines from the combined SNela 
and CMBR shift parameter data. In the gray area the scal- 
ing solution is unstable regardless of a. If the scalar field 
is tracking already in the nucleosynthesis epoch, there is a 
tension with the amount of early quintessence and the nucle- 
osynthesis limit. Conservatively, this limit would translate 
to A > 6.3 [4]. 

lates the gradient of the gravitational potential to 
the perturbations in matter density. The result in 
our case is 

G. _ -x' + + + x^[2{l + e)ip-l) + yj 
G x2[4 + ^I{5^i - 8)] - /i2[6(l + e)ifi - 1) + yV 

consistently with Ref. [27] . Though not obvious from 
the formula, it equals one when the coupling goes 
to zero. Generally, when /i is of order one, then one 
expects the G'^/G* to be of roughly of order one 
as well. As claimed in Ref. [27], one has to assume 
" an accidental cancellation" to satisfy the bound for 
G'^/G* in the presence of significant GB contribu- 
tion to the energy density. This means that we 
have to fine-tune the coupling parameter a in or- 
der to eliminate the time variation of the effective 
gravitational coupling, typically with the accuracy 
of 0.01. Then the Newtonian limit in general ex- 
hibits time- varying G* , but at present days this G* 
appears to us as a constant. We have reported the 
cosmological results in such a case in Table 4. The 
best-fit values of per effective degree of freedom, 
Xdof' slightly better than in the ACDM case. 



* Such a tight bound might not be so problematical if one 
takes into account that cosmological variations of G and 
other gauge-couplings might be different from the ones we 
measure on Earth or within our Solar system [29]. 



4 



Data set 


ACDM 


RQg4> model 




^dof 








A 


a 


SNela 


1.142 


0.314 


1.146 


0.42 


5.1 


32.3 


SNeIa+7^ 


1.144 


0.277 


1.141 


0.44 


5.2 


33.8 



Table 1 

Best-fit values for ACDM model compared with fits of the 
coupled scalar field model for some parameter values when 
the coupling a is set in order to fix the present time variation 
of G* to zero. The degrees of freedom in the first row are 
157 — d, and in the second 158 — d, where d = 1 for the 
ACDM and d = 3 for the R%g4> models. 

Interestingly, G* determines the cosmological 
evolution of matter inhomogeneities at subhorizon 
scales. The matter overdensity 6 evolves according 
to [20] 



2HS ^ AirG^pS. 



(8) 



From this equation we can make the important con- 
clusion that the evolution of matter inhomogeneities 
at small scales is independent of wavelenght. Thus 
we expect, given the same primordial spectrum of 
perturbations, the main difference in the matter 
power spectrum (when compared to ACDM) to be 
in the normalization. Note that to arrive at this 
result, we have not assumed specific forms for the 
coupling or the potential, but only that the scalar 
appearing in the action (2) is not very massive. Our 
numerical solutions of the whole perturbation sys- 
tem confirm that the approximation Eq.(8) indeed 
is good at subhorizon scales. 

The effective gravitational constant G» can in 
principle diverge. We find that in the example 
model presented here, G» typically diverges in the 
future, and for low matter densities this can happen 
even before a = 1. It is unclear what happens at 
such point, since the linear approximation certainly 
breaks down near the (what would be) singular- 
ity. Matter perturbations will at least for a while 
grow explosively. It is possible that consequently 
the de Sitter phase will not be then reached, which 
could help to define the S-matrix in string theory. 
This divergence can be related to previous stability 
considerations of these models (for recent studies, 
see Refs. [28,30]). This far stability conditions for 
these models have been derived only in the case 
that the field perturbations decouple from other 
fiuids. However, in dark energy cosmologies such as 
here one should take into account both matter and 
the corrections to Einstein gravity. We provide a 
clear indication that the scalar action in its vacuum 
form determines the stability also in the case that 



Pm 7^ 0. The action for the potential <f> in vacuum 
features an effective propagation speed [22] 

_ _a;2[4 + ^(5^ - 8)] -t- p^[6{l + e){fi -1) + y] 
(m-1)[3m2-4(/.-1)x2] ^> 

from which we can see that the linearized matter 
perturbations diverge exactly at the points where 
this propagation changes its sign [3- 

To summarize, we studied a string-inspired low 
energy action where a coupling between a scalar 
field and the GB invariant is present. Solving the 
full system of equations which describe both the ho- 
mogeneous universe as well as its first order pertur- 
bations we found several new insights to GB dark 
energy cosmologies in general, which could be rele- 
vant to future investigations of viable string-theory 
low energy effective actions. We showed that the 
background universe can present an attractor solu- 
tion which features late transition from a scaling 
era to acceleration triggered by the GB coupling. 
This mechanism might be seen to alleviate the co- 
incidence problem. In addition, the model can con- 
sistently explain a presently ongoing but transient 
phantom era. This background expansion can ex- 
hibit a good a accordance with cosmological and as- 
trophysical data. The evolution of matter perturba- 
tions is scale- invariant at small scales in the presence 
of the GB term, and thus the shape of the matter 
power spectrum is retained. Hence, the latest data 
from the CMBR anisotropics as well as LSS can 
show agreement with these models. However, using 
a combined set of present days cosmological obser- 
vations it is possible to constrain the parameters of 
the theory tightly. 
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